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ABSTRACT 


Drake has suggested that apparent multimillisecond 
periodicities in pulsar signals may correspond to vibrations 
of neutron stars. We have reexamined stellar models in the 
appropriate density range with an improved equation of state, 
and we conlude that suitable stable structures are unlikely. 
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At a recent seminar, Drake (Drake 1968) haB reported a 

superposition of 10 and 15 millisecond periods on the pulses 

received from the pulsars AP 2015 and CP 1919. Ha ia reported 

to have suggested that these periods represent the intrinsic 

vibrational period of a neutron star. However, neutron star 

periods much in excess of one millisecond would be very unlikely 

in a neutron star (Tsuruta, Wright, and Cameron 1965) . A period 

of order 10 milliseconds would be roughly characteristic of a 

13 , 3 

body with a mean density near 10 gm/cm , two orders of magnitude 

smaller than the mean density of a typical neutron star. The 

composite white dwarf-neutron star model curve of Tsuruta and 

13 

Cameron (1966) does show a subsidiary mass peak between 10 
and 10*"^ gm/cm'* . However, because the treatment of the physics 

m 

in this density range was extremely crude, no significance was 
claimed for this peak by the above authors. However, the 
suggestion by Drake raises an important question as to the 
possible existence of stable structures in this density range. 

In this note we report on an investigation of this question with 
an improved but still crude treatment of the relevant physics. 

The following physical phenomena play a key role in the 
equation of state. In the density range up to about 3 x 10 11 gm/cm 3 
only ions and electrons are presert under cold degenerate conditions.' 
In this range the electron Fermi level approaches 30 MeV, and 
electron capture reduces the charge to mass ratio of the nuclei. 
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The neutron binding energy of the nuclei approaches zero; the 
point where it becomes zero has been denoted the "neutron drip 
line" by J.A. Wheeler. At still higher densities and electron 
Fermi energies, electron capture continues, and neutrons leave 
the nuclei. Gradually a degenerate neutron gas is established 
in the volume between the nuclei, and a significant neutron 
Fermi energy is established. At this stage the nuclei will 
contain unbound neutrons, since neutrons 'can only be removed from 
the nuclei with energies in excess of the neutron Fermi level. 

Of course, neutrons will be readily exchanged between the free 
neutron sea and the nuclear interior, but the density of unbound 
neutrons in the nuclear interior will considerably exceed that 
outside, since the Fermi level of the neutrons inside the nucleus 
will be much greater than that ootside, and thus the number of 
momentum states per unit energy interval near the Fermi surface 
will be much greater inside than outside. 

We have treated the physics in the following very approximate 
way. We write tyft=M-A for the mass excess of a nucleus. We use 
the simplest possible form of the mass equation: 

^f)(A,Z) Z+ < ^ 1 (A-Z)+OA+0(A-2Z) 2 /A+yA 2 ^ 3 +6Z 2 /A 1/3 (1) 

Here and ^5^ are the mass excesses of the proton and neutron, 

and the last four terms on the right hand side represent respectively 
the volume binding energy, the volume symmetry energy, the surface 
energy, and the coulomb energy. This equation is applicable to 
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nuclei with sharp edges. In what follows we shall make some 
approximations valid for large mass numbers A. 

The beta decay energy is 

Q -97HZ,A) - 0»7(Z+1 # A) (2) 

This leads to the expression 

Z -a -( %,-%,) + 46 , - 4?/A - 6/a 1 / 3 

A 83+ 26A 2/3 

We use values of the coefficients given by Green (1954) : 
a » -15.756 MeV, 0 - 23.694 MeV, Y * 17.794 MeV, and 
6 = 0.7103 MeV. It may be seen that the last two terras of the 
numerator of (3) are small compared to the preceeding ones, 
and that the second term in the denominator is fairly small 
compared to the first. Hence for a given value of the beta 
decay energy the ratio Z/A is not strongly dependent on mass 
number. Taking as typical A'vlOO, we have 


A 220 

The neutron binding energy is 

B = r /0 n + < %(Z,A-1) -^1(Z,A) 


This leads to the expression 


4/3 

Z 3 (B+a+fl) + 26/A 
A2* 123 


?■> f-2.6/? 

+ 6 * 2/3 


(4) 


(5) 


( 6 ) 


It may be seen that the last terms in the denominator and 
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numerator of (6) are small compared to the preceeding terms, and 
hence for a given neutron binding energy the ratio Z/A is not 
strongly dependent on mass number. Again taking as typical 
A 100 , we have 


Z 2 . B+7.9524 
A J 100 


(7) 


The neutron drip line (B=0> occurs at Z/A = 0.282, 

according to (7) . From (4) , the corresponding beta decay energy 

(or, equivalently, electron capture threshold) is Q = 31.96 MeV. 

An equation of state can now be constructed, with the pressure 

and density parametrically depending on the electron Fermi energy 
P 

E . The latter quantity specifies n and P , the electron number 
density and pressure. The total number of protons (inside 
nuclei) per cm^ is equal to n , for charge equality. We rewrite 
(4) as 


Z = 94 - El (8) 

A 220 


This specifies the number of neutrons in nuclei per cm 3 , and 

hence the mass density of the nuclei. The ionic pressure can be 

F 

neglected. If E^ < 31.96 MeV, the equation of state is thus 

F 

completely specified. If E e > 31.96 MeV, then free neutrons 

F 

are also present. We calculate the neutron Fermi energy E 

n 

from a modification of (7) : 


Z 2 -E„ + 7.9524 

— T ■ ~ 1 ** 

A 2 100 


(9) 
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From E n we also calculate the number density of free neutrons in 

the internudear region and the neutron pressure P n , which will 

also apply to the region occupied by nuclei. We assign the 

nuclei a density 4 X 10 14 gm/cm 3 , calculate the fractional 

volume occupied by the nuclei, and hence find the mass of the 

neutrons in the internudear region. The total pressure P ■ 

P e + P n . The equation of state is thus completely specified. 

It may be noted from equation (8) that Z/A — 0 as - 94 

e 

MeV. This indicates, since the total number of protons per 

F 

unit volume is fixed by the value of E e , that ever increasing 
numbers of neutrons are absorbed into the nuclei. If individual 
nuclear charges become too small, such nuclei will fuse by 
pycnonudear reactions. Thus the fractional volume occupied 
by the nuclei approaches unity as the total density approaches 
4 X 10 14 gm/ cm3 and as E^ - 7.9524 MeV. 

This equation of state is obviously not a good approximation 
near a density of 4 X 10^ 4 gm/cm 3 and certainly cannot be used 
beyond this density. However, this limit is quite sufficient 
for the present investigation, which is intended only to bring 
out the essential physical behavior of the equation of state. 

A more refined investigation should take into account higher terms 

i 

in the nuclear symmetry energy, the rounded edge of the nucleus, 
corrections to the surface energy when a neutron sea is present , 
corrections co the coulomb energy due to the electrons, formation 
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of a crystalline lattice, and nucelar potential energies of the 
free neutrons. 

With this equation of state the equations of hydrostatic 
equilibrium were integrated to obtain stellar models. Figure 1 
shows the relation between stellar mass and central density, 
showing the well known dwarf mass peak and the region up to 
a central density of 4 X 10^ gm/cm^. It may be seen 
that there is no second mass peak in the "higher density region, 
although the mass starts rising slowly above log p=13.5. 

In the Tsuruta-Cameron procedure the neutronization of 
nuclei was assumed to create new nonrelativistic particles 
as the matter was compressed, thus assuring a large value of 
the adiabatic index P^. Their secondary peak resulted from 
these assumptions. In this improved procedure, the neut- 
ronization process initially creates many nonrelativistic 
neutrons, but these are progressively absorbed on nuclei as 
compression occurs at the higher densities. Hence remains 
low The values of fj* obtained for the equation of state de- 
scribed here are shown in Figure 2. It may be seen that 
attains a maximum value of 1.24 in the neutronization regio.i at 
log p-“ 12.8. To create a secondary meass peak, major portions 
of a star would require ^ > 4/3, which is not possible here. 
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FIGURE CAPTIONS 


Figure 1. Mass of a degenerate star plotted as a function of 

central density, in the range below nuclear density. 

Figure 2. Adiabatic Index in the density range up to 

nuclear density corresponding to the equation of 
state developed in this paper. 













